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SUMMARY 


A  theory  is  developed  for  the  one-dimensional  motion  of  a  cylindrical 
gas  bag  used  as  an  impact  cushion.  The  effect  of  shock  waves  in  the  gas 
as  well  as  stress  in  the  bag  skin  is  considered.  The  applicability  of  the 
theory  to  landings  both  in  an  atmosphere  and  on  the  moon  is  discussed  and 
the  regime  of  validity  of  the  theory  is  presented.  The  use  of  a  series 
expansion  for  computing  shock-wave  properties  in  the  analysis,  the  strong 
shock  approximation,  and  the  exact  shock  relations  are  compared  and  dis¬ 
cussed.  The  regime  of  physical  parameters  for  which  both  the  wave  model 
and  the  series  expansion  are  valid  is  presented.  The  method  of  application 
of  the  theory  to  impact  problems  is  outlined. 


INTRODUCTION 


Over  the  past  five  years  there  has  been  f  growing  interest  in  impact 
attenuation.  Considerable  experimental  and  sor.fi  theoretical  work  has  been 
done  on  the  study  of  impact.  Drop  tests  have  teen  performed  at  the 
University  of  Texas  (refs.  1  to  5)  using  numerous  impact  attenuation 
techniques,  including  crushable  structures,  foamed  plastics,  and  inflated 
cylindrical  bags.  Most  of  that  work  was  performed  at  low  impact  speeds. 

The  use  of  inflated  bags  for  cushioning  the  impact  of  ejected  aircraft 
pilot  compartments  has  been  investigated,  as  well  as  impact  attenuation 
associated  with  atmosphere  entry  vehicles  (refs.  6  to  8).  A  theory  of  low- 
speed  impact  of  various  shaped  gas  bags  has  been  developed  by  Esgar  and 
Morgan  (ref.  9)*  Martin  and  Howe  have  developed  a  uniform  compression 
theory  for  the  impact  of  inflated  spheres  (ref.  10)  as  well  as  a  theory 
which  accounts,  in  an  approximate  way,  for  wave  motion  in  the  inflating 
gas  (ref.  ll) .  In  these  last  two  papers,  a  two-dimensional  unsteady 
problem  was  simplified  in  various  ways  to  achieve  a  workable  analysis. 

The  problem  in  the  present  analysis  has  a  simple  one-dimensional 
geometry.  For  this  reason,  fewer  assumptions  need  be  made  in  the  develop¬ 
ment  of  the  theory.  Basically,  the  problem  has  one  space  dimension  and 
one  time  dimension  as  independent  variables.  Ihe  analysis  is  valid  up  to 
the  time  that  the  payload  comes  to  rest.  In  solving  this  problem,  it  is 
necessary  to  couple  the  effects  of  the  unsteady  wave  motion  in  the 
inflating  gas  with  those  of  the  structural  behavior  of  the  bag  skin. 


The  latter  aspect  of  the  problem  could  be  exceedingly  difficult  in  itself 
if  one  attempted  to  solve  it  as  an  unsteady  elasticity  problem  with  a 
moving  discontinuous  load.  This  would  introduce  two  space  and  one  time 
dimension,  which  when  combined  with  the  unsteady  gas  dynamics  might  make 
solving  the  problem  a  virtually  hopeless  task.  The  waves  in  the  skin 
would  be  coupled  with  the  waves  in  the  gas  although  the  two  would  not 
move  at  the  same  speed  and  the  motion  of  the  payload  would  be  exceedingly 
complicated. 

Instead,  we  adopt  a  simple  model  for  the  bag,  considering  it  to  be 
a  cylindrical  membrane  with  high  sound  speed,  and  devote  most  of  our 
attention  to  the  wave  motion  in  the  inflating  gas.  The  interaction  of 
families  of  gas  waves  is  examined,  and  a  reasonable  gas  wave  model  is 
chosen  for  analysis.  The  conditions  under  which  the  model  is  valid  are 
examined  for  impact  on  a  planet  with  an  atmosphere  as  well  as  on  an  air¬ 
less  surface.  The  property  relationships  across  the  shock  waves  in  the 
gas  are  expressed  in  three  ways:  by  the  series  expansion  method  and  the 
strong  shock  approximation  which  have  the  advantage  of  analytical  sim¬ 
plicity,  and  the  exact  method.  The  conditions  that  limit  the  use  of  the 
first  two  methods  are  discussed.  Relationships  for  maximum  acceleration, 
stopping  distance  (stroke),  bag  stresses,  and  pressures  are  derived  from 
the  theory. 


SYMBOLS 


A  cross-sectional  area  of  cylindrical  bag,  and  coefficients  in 

equation  (48) 

c  sound  speed  in  the  inflating  gas 

c  dimensionless  sound  speed  defined  by  equation  (8) 

F  forces  in  sketches  (a)  and  (b) 

g  gravitational  constant  at  the  earth's  surface 

K  defined  by  equation  (ll) 

m  mass 

mass  of  sides  of  bag 

n  allowable  number  of  g's  acceleration 

N  force  per  unit  length 

p  "  pressure 
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P 


cross-sectional  radius  of  cylindrical  bag 

gas  constant 

time 

dimensionless  time  defined  by  equation  ( l4) 
temperature 

velocity  in  positive  y  direction 
dimensionless  velocity  defined  by  equation  (8) 
instantaneous  volume  of  gas  bag 
shock  speed  relative  to  gas  at  state  i 
defined  by  equation  (4-9) 

distance  normal  to  impact  surface  (positive  outward) 
dimensionless  distance  defined  by  equation  (8) 

defined  by  equation  (23) 

ratio  of  specific  heats  of  inflating  gas 
mass  per  unit  surface  area 
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Subscripts 


o,  1,1  states  shown  on  sketch  (d),  and  subscripts  on  coefficients  in 
a,  3  J  equation  (48) 


Ol,  12,1 
23  J 

a 

b 


d,  e,  f 


g 


shock  connecting  states  0  and  1,  1  and  2,  2  and  3 

point  on  sketches  (c)  and  (d)  or  atmosphere,  axial 

point  on  sketches  (c)  and  (d)  or  side  walls  of  bag 

point  on  sketches  (c)  and  (d)  or  circumferential 

points  on  sketches  (c)  and  (d) 
inflating  gas 


i 


condition  ahead  of  a  shock 


i 


J  conditions  behind  a  shock 

p  payload 

bottom  bottom  of  bag 


ANALYSIS 

The  Equation  of  Motion 


Consider  a  nonslender1  cylindrical  bag  of  height  y0  and  cross- 
sectional  area  A.  The  bag  is  inflated  to  a  pressure  pQ.  On  top  of  the 
bag  is  a  payload  of  mass  mp.  The  pressure  on  the  outside  of  the  payload 
is  pa.  It  is  assumed  that  the  bag  is  stiff  enough  to  remain  cylindrical 
under  p0  -  pa-  The  system  moves  along  the  axis  of  the  right  cylinder 
with  speed  | Uq  |  and  strikes  a  flat  surface  normal  to  the  cylinder  axis. 
During  the  impact,  the  pressure  in  the  inflating  gas  at  the  payload  is  p 
(at  height  y  above  the  surface).  A  free  body  diagram  of  the  payload  is 
shown  in  sketch  (a)  (where  it  is  assumed  that  the  bottom  of  the  payload 
is  flat  and  always  normal  to  the  axis).  The  momentum  theorem  applied  to 


Po 


Sketch  (a) 

that  free  body  is,  with  gravity  forces  neglected, 

mpY  ■  (P  -  Pa)A  “  Fl  (!) 


XA  nonslender  bag  is  defined  as  one  not  subject  to  long  column  type 
buckling.  For  the  case  of  a  circular  cylindrical  bag  and  low  impact  speed, 
it  appears  that  if  the  diameter  is  less  than  about  75  percent  of  the 
initial  bag  height,  the  gas  bag  buckles  and  is  unstable  (ref.  5). 
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A  free  body  diagram  of  the  sides  of  the  bag,  excluding  the  gas,  is  shown 
in  sketch  (b).  In  time  dt,  the  skin  changes  momentum  by  the  amount  mSvdu 


as  a  result  of  a  change  in  velocity  (where  msv  is  the  bag  mass  in 
motion).  In  the  same  time,  the  momentum  change  due  to  a  change  in  mass  is 
dm6v(u) .  The  total  momentum  change  is  m8Vdu  +  dmsv(u)  »  (Px  +  F2)dt, 
where  shear  forces  due  to  gas  or  atmosphere  are  neglected  or 


~  (mBVu)  -  FX  +  F2  (2) 

The  forcr  F2  brings  the  mass  dmBV  to  rest  in  time  dt.  Here  the 
unrestrictive  assumption  has  been  made  that  the  tension  goes  to  zero 
(the  buckling  stress  of  the  bag)  at  the  top  of  the  element  dmsv. 


Fadt  *  u  dmgv  (3) 


Therefore  from  equations  (2)  and  (3) 

Fi  =  msv  ~  (4) 

From  equations  (l)  and  (4) 

#  (mp  +  msv )  u  (p  -  Pa)A 


(5) 
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If  is  the  mass  o'  the  sides  of  the  hag,  then 


then 


®sv  -  ~ 


(“P +  *b  57)  «  t  ■ (p  ■  P*)A 
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If  K  is  defined  as 


(mp +  |f  a  poA(p  -  pa) 


y0P0A  P0V0  _  mg  KTq  _  mg  c0£ 
mpUo2  np\ioa  nip  Uq2  mp  7 


then 


By  definition 


dg_  K(F  -  pa) 

<3y  u[l  +  (n^y/nip)! 


dt  _^L 
dy  u 


t -±aL 

y0 


where 
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Now  to  integrate  equations  (12)  and  (13) ,  we  need  initial  conditions  and 
a  description  of  p  and  pa.  For  very  low  inqpact  speeds  (high  cT0)  it  is 
probably  adequate  to  describe  p  by  the  isentropic  relation 

pV7  -  1  (15) 

where 

=  7  (16) 

'0  J'o 

However,  for  high  impact  speeds,  equation  (15)  is  not  adequate  because 
shock  and  rarefaction  waves  move  through  the  gas,  and  the  pressure  must 
be  calculated  accordingly.  This  will  be  considered  in  the  next  section. 

Similarly,  we  need  to  express  pa.  On  the  moon,  it  is  zero.  In  an 
atmosphere,  the  entire  payload  and  inflated  bag  combination  moving  toward 
the  surface  has  a  flow  field  around  it,  a  wake,  and  a  stagnation  region. 
When  the  bag  strikes  the  surface,  the  flow  is  altered  (including  the  wake) 
and  the  pressure  in  the  immediate  neighborhood  of  the  system  varies  with 
space  and  time.  A  detailed  treatment  of  this  external  flow  field  is 
obviously  very  complicated.  Although  various  approximations  to  the 
external  flow  could  be  adopted,  we  adopt  the  simplest  one  possible  and 
say  that  the  external  pressure  on  the  payload  is  represented  by  an 
effective  constant  average,  pa.  We  now  move  on  to  the  calculation  of  ^ 
by  considering  the  wave  motion  in  the  inflating  gas. 


The  Wave  Model  in  the  Inflating  Gas 

Sketch  (c)  represents  the  most 
general  wave  pattern  in  the  inflating 
gas.  The  curve  abfd  is  a  trace  of 
the  payload  motion  in  the  space -time 
plane.  The  heavy  curves  represent 
the  traces  of  shock  waves  and  the 
dashed  lines ,  those  of  rarefaction 
waves.  At  time  zero,  the  bottom  of 
the  bag  strikes  the  ground  and  shock 
wave  ob  propagates  upward  (analogous 
to  the  impulsively  started  piston). 

At  or  near  time  zero,  stress  relief 
in  the  bag  skin  changes  the  motion  of 
the  payload,  ab,  causing  a  rarefaction 
to  propagate  downward  from  ab.  That 
rarefaction  intersects  the  initial 
upward  traveling  shock  at  e.  The 
shock  is  modified  by  the  rarefaction 
and  its  trajectory  eb  is  curved. 


0  c 


t 

Sketch  (c ) 


;  The  shock  then  reflects  from  the  "bottom  of  the  payload  at  h  and  the 

f  trajectory  of  the  reflected  shock  he  is  also  curved.  The  reflected 

1  head  of  the  initial  rarefaction  wave  crosses  the  reflected  shock  and 

\  arrives  at  the  payload  at  f .  In  general,  the  computation  of  the  pressure 

at  the  payload  depends  on  what  happens  in  the  entire  wave  system  and  is  a 
i  difficult  computation  for  the  complete  wave  system  shown  in  sketch  (c). 

On  the  other  hand,  if  the  change  in  forces  (and  motion)  on  the 
payload  due  to  stress  waves  in  the  skin  is  negligible  compared  with  that 
due  to  pressure  waves  in  the  gas,  the  wave  pattern  is  much  simpler 
(sketch  (d)).  Uiis  condition  is  approached  if  pg  -  pa  »  pQ  -  pa  and 

is  exact  if  pQ  =>  p&  (no  overpressure). 

Even  if  pa  is  zero  ( lunar  impact )  the 
approximation  is  valid  provided  p2 
is  much  larger  than  pQ  (corresponding 
to  small  c0).  We  will  adopt  the 
simple  wave  pattern  of  sketch  (d)  and 
focus  our  attention  only  on  cases  where 
p2  -  pa  >  10(1  “  Pa)i  the  arbitrary- 
factor  10  is  considered  conservative . 

Thus  the  paths  ab  and  ob  are  straight 
lines.  Paths  be  and  cd  may  be 
curved,  but  they  do  not  influence  the 
path  bd  or  the  computation  of  the  * 

pressure  along  bd  except  in  locating 
point  d.  Subsequently,  be  and  cd 
will  be  treated  as  straight  lines  and 
will  be  discussed  in  detail. 

If  the  payload  comes  to  rest 

Sketch  (d)  before  the  point  d  is  reached,  the 

entire  pressure  computation  is  that 
of  a  simple  wave  (ref.  12,  p.  4l3)  .  In  the  simple  wave  region 
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and  from  isentropic  relationships  (noting  that  u®  =*  -1  since  the  slope 
of  abd  is  continuous  at  t>) 

2 7  27 


JL 
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Substituting  equation  (19)  into  (12)  yields 


(19) 


du  _ 
dy  " 


KP£ 


27 

jy-i  _  1 

"  pa 


Ki+=?y, 

For  convenience  equation.  (13)  is  repeated 


dt  =  1 
dy  u 

The  boundary  conditions  on  equations  (20)  and  (13)  are  at 


(20) 


(13) 


at  y  -  yb 


u  =  -1 
t  =  tb 


(21) 


Equations  (20) ,  (l3)>  and  (2l)  comprise  the  set  of  differential 
equations  and  their  'boundary  conditions  whose  solution  is  the  motion  of 
the  payload  during  impact.  Before  describing  ways  to  obtain  p2,  ~c2,  y^, 
and  it  is  instructive  to  integrate  equation  (20)  formally  for  the 
case  =>  0. 


Payload  Velocity  During  Impact  Without  Atmosphere 


Variables  in  equation  (20)  are  at  once  separable  (note  that  pa  is 
set  to  zero) 


dy 


u  du 


1  +  ^7 

J 


(22) 


mp 


10 


Denoting 


a 


2c2  +  1-7 
25=, 


3  a  -fr  -  D 

^  r» 


250 


T  B 


mp 


> 


(23) 
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and  noting  that  p2  and  c2  are  constants  determined  by  initial  conditions 
(as  will  be  shown  subsequently)  we  integrate  equation  (22)  subject  to 
boundary  conditions  (21).  This  gives  (since  n  ^  -1  or  -2) 


y(u)  -  (i  ♦ 


[  T 

F(a  +  -  1 

+  2 

V-  +  l 


1 

T 


(24) 

Equation* (24)  represents  the  motion  of  the  payload  during  intact  in  a 
vacuum.  It  does  not_appear  to  be  readily  integrable  again  in  combination 
with  (13)  to  yield  t(3r) .  Note  that  at  u  =  0, 
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tt 


yinin  "  1 7  +  ybJexPi 


/  T 

r  p+2  , 

a  -  1 

Wp2 

\i  +  2 

n  +  1 


1 

T 


(25) 


For  large  'c0,  ”2  is  large  and  a  -*  1,  so  that 

ymin  yb 


(26) 


In  the  special  s£lution  of  .equation  (25)  and  in  the  general  problem, 

P2>  yb#  6111(1  ^b  are  still  not  determined.  These  will  now  be  deter¬ 
mined  from  normal  shock  relationships. 
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Normal  Shock  Relationships 


Having  adopted,  a  simple  wave  model  and  included  it  in  the  differen¬ 
tial  equation  of  motion  of  the  payload,  we  now  compute  initial  conditions 
Just  inside  the  simple  wave  region  at  b  (sketch  (d))  by  three  methods. 

The  first  method  is  a  series  expansion.  It  has  the  advantage  of  simplic¬ 
ity,  and  leads  to  explicit  physical  relationships.  However,  its  validity 
is  limited  to  certain  regimes  of  impact  conditions.  The  second  method 
involves  the  use  of  exact  shock  relationships,  and  the  third  is  the 
strong  shock  approximation. 

Series  expansion.-  The  series  expansion  relationships  for  properties 
across  a  moving  shock  are  valid  as  long  as  the  shock  is  not  too  strong. 

In  the  present  paper  the  restriction  on  shock  strength  corresponds  to  a 
restriction  on  c‘0.  It  will  be  shown  subsequently  that  if  ~0  >  1,  there 
will  be  an  error  in  the  series  expansion  computation  of  p2  and*”^,  of 
10  percent  or  less  (for  y  =  1.4l) .  This  corresponds  to  a  maximum  value 
of  p2  of  roughly  9  for  use  of  the  series  expansion. 

In  addition  to  this  restriction  associated  with  the  series  expansion, 
there  is  the  further  restriction  associated  with  the  wave  model  adopted 
in  sketch  (d) .  For  the  wave  model  to  be  adopted  it  was  required  that 

P2  “  ^a  >>  ?o  “  Pa  =  1  "  ^a 

If  pa  =  0,  and  since  the  upper  limit  on  p2  a  9,  it  is  seen  that  the 
inequality  is  not  really  satisfied  in  the  strong  sense  by  solutions  using 
the  series  expansion.  Thus  the  series  expansion  should  not  be  used  for 
impact  without  an  atmosphere;  the  exact  shock  relations  or  strong  shock 
relations  should  be  used  instead.  On  the  other  hand,  if  pa  =  1,  the 
inequality  is  always  satisfied,  the  wave  model  in  sketch  (.d)  is  applicable, 
and  the  series  expansion  can  be  used  as  long  as  c0  is  greater  them  unity 
as  mentioned  above  (but  less  than  an  upper  bound  that  will  be  shown 
sub  sequent  ly) . 

From  sketch  (d),  it  is  seen  that 

Uq  =  -1 
Ux  =  0 

u2  =  -1 
u3  =  0 


(27) 
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The  polynomial  expansions  of  shock  relations  in  terms  of  shock 
strength  are  (ref.  12,  p.  1011) . 


ci  -  c0*+  7-  1  (ui  -  Uq)  +  .  .  .  [  0]  (ui  -  Uq)3  a  cc 


7  -  1 


(28) 


C2  =  Cl  -  ~~  (ua  -  Ui)  +  .  .  .  [0](u2  -  Ui)3  *  Co  +  7  -  1  (29) 


w, 


-  up  +  Cp  +  ui  +  cj.  _  jjgo  +  y  -  3 


shock  oh 


Trajectory  ob  is  (integrating  eq.  (30)) 


(30) 


A 

5 

7 

k 


_  _  ^4c  +  r  -  3^ 


(3D 


and  of  course  the  trajectory  of  the  payload  ab  is 

y  =  1  -  t 

From  equations  (31)  and  (32) ,  solve  for  y  =  y^  and  t  =  tt, 


(32) 


_  _  ^cp  +7-3 
=  4co  +  7  +  1 


(33) 


\ 2 — 
D  4C0  +  y  + 


(3*0 


It  is  possible  to  stop  at  this  point  and  compute  p2  which  would 
give  all  the  information  needed  to  start  the  integration  of  equations  (20) 
and  (13).  However,  we  also  want  to  have  some  understanding  of  the  extent 
of  the  simple  wave  region  so  that  we  know  whether  or  not  the  shock  cd 
reflects  off  the  payload  before  the  payload  comes  to  rest.  If  that  were 
the  case,  we  would  be  out  of  the  single  wave  region,  and  the  simple  model 
would  be  invalid  for  that  part  of  the  motion  after  the  reflection  at  d. 

An  approximate  knowledge  of  the  location  of  point  d  will  satisfy  pre-  • 

sent  purposes;  so  in  the  interest  of  simplicity,  straight  lines  will  be 
used  for  both  be  and  cd.  The  computation  of  the  pressure  and  motion 

» 


in  the  simple  wave  along  hd  is  not  affected.  Only  the  location  of  d 
is  affected.2  The  straight  line  he  is  determined  by  integrating  the 
equation 


TTj.  -  “j.  +  Ue  -  "c2  kTS0  +  37-1 
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(35) 


which  leads  to 


? .  .(m-> ♦ 


(36) 


When  y  3  0,  t  3  tc 


t„  = 


16(2cq  +7-1) 


(4c0  +  7  +  l) ( 4c0  +37-1) 


State  3  &t  the  ground  is  given  by 


c3  =  72  +  7—g’1-  (U3  -  us)  =  c0  +  ^  (7  -  1) 


(37) 


The  slope  at  the  ground  of  the  reflected  shock  cd  is 


/d£\  _  Ua  +  c2  +  U3  +  c3  =  4cQ  +  57  -  7 

,dV23  ^ 


(38) 


The  equation  of  the  straight  line  cd  is  {integrating  (38)  and  saying 
3  0  at  ^  ~  ^q) 


t  =  4  _  +  l6(gcQ  +  7-1) 

+  57-7  y  (4c0  +  7  +  l)(4c’0  +  37-1) 


(39) 


2It  can  be  shown  that  be  curves  right,  but  that  cd  may  curve  left 
or  right.  If  both  be  and  cd  curve  right,  the  straight  lines  give  a 
conservative  location  for  d  inasmuch  as  solutions  will  be  obtained  for 
which  the  payload  comes  nearly  to  rest  before  the  approximate  d  is 
reached  in  advance  of  the  exact  d.  Also,  if  be  curves  right  and  cd 
curves  left,  one  counteracts  the  other,  and  the  straight  line  is  probably 
satisfactory. 


The  intersection  of  (39)vlth  the  payload  trajectory  obtained  by  the  « 

simultaneous  integration  of  equations  (20)  and  (13)  gives  point  d. 

Next  p2  is  evaluated  by  the  series  expansion  (ref.  12,  p.  1011) 


Ill 

,  .  rtai  -  it) 
*0 

,  7(7 
•  + - 1 

+  1) 
i 

<N?)  *••• 

(40) 

£1- 

4702  +  4^07  + 

0 

72  +  7 

(41) 

Po 

4^o 

w 

II 

*  -  <*  ;,•*) 

)  2±n 

+  1) 

T~ 

(42) 

Pa  = 

4c02  +  4(27  - 

l)c0  + 

4y2 

-37+1 

(43) 

Pi 

4c02  +  4(7 

-  l)c0 

+  (7 

-  I)2 

From  (h-l)  and  (43)  since  pD  =  1 

_  ^  [4c02  +  4(2y  -  l)c0  +  by2  -  3y  +  l][4c02  +  4yc0  +  7(7  +  l) ] 
2  4c02[4c02  +  4(y  -  l)c0  +  (7  -  l)2] 


Thus  the  problem  is  completely  formulated  by  the  series  expansion  method. 

To  solve  a  given  case,  one  specifies  ^0,  7,  K,  m^/mp,  and  p&.  Then 
equations  (20)  and  (13)  are  integrated  simultaneously  subject  to  boundary 
conditions  (21),  (33)  ,  and  (34)  where  ‘c2  and  p2  are  calculated  from 
equations  (29)  and  (44) .  Each  example  should  be  checked  to  be  sure  that 
u  =  0  by  the  time  the  intersection  d  has  been  reached  (by  use  of 
eq.  (39))  • 

Exact  shock  relations.-  If  “q  <  1,  it  was  stated  (for  7  =  1.4l) 
the  exact  shock  relations  or  the  strong  shock  approximation  must  be  used 
for  the  computation  of  p2  to  be  accurate.  It  was  also  mentioned  that 
for  impact  without  an  atmosphere,  p2  must  be  10  or  more  for  the  wave 
model  to  be  valid.  This  requires  that  c’0  be  less  than  1  and  thus  the 
exact  shock  relations  need  be  used  to  determine  the  starting  conditions 
(c2>  p2,  y^,  ^j)  for  integration  of  the  differential  equations  (20)  t 

and  (13).  In  addition,  the  assumed  straight  line  trajectories,  be  and  cd, 
should  be  computed  by  the  exact  shock  equations. 


<  UA  f--4- 


•  > 
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In  sketch  (e) ,  the  shock  represented  by  the 
solid  line  is  moving  into  the  gas  at  state  i  with 
speed  ±W£ j  relative  to  the  gas  at  state  i.  The 
gas  on  the  backside  of  the  shock  is  at  state  j . 

The  plus  sign  on  denotes  upward  traveling 

shocks  (with  respect  to  gas  at  state  i)  and  the 
minus  sign  denotes  downward  traveling  shock  waves. 

The  sub scripts  i  and  j  used  to  denote  states  i 
and  J  will  eventually  take  on  the  pairs  of  val¬ 
ues  01,  12,  23  when  the  shock  relations  are  applied 
to  the  shocks  connecting  states  0  and  1,  1  and  2,  and  2  and  3  in 
sketch  (d) .  The  exact  shock  relations  made  dimensionless  by  equations  (8) 
and  (l4)  are  (see  ref.  12,  pp.  1001,  1002) 


Sketch  (e) 


Solving  equation  (45)  for  Wij/ci  and  combining  the  results  with 
equation  (47)  yields 


where 


and 


» 


t 


xj4  +  A3X13  +  A2X12  +  AiXi  +  A0  =  0 


*i 


Aq  = 


Ai 


A2 


0 


'(y  +  D2 

.  47Ci2 


A3 


(y  +  D2 

2y*i 


j 


(48) 

(49) 


(50) 
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The  application  of  the  equations  to  the  shocks  connecting  the  ij  '  ( 

states  01,  is  as  follows.  It  is  noted  that  Uq  =  -1,  and  uj.  =  0. 

Equation  (48)  is  solved  for  xQ  hy  use  of  equation  (50)  where  7:1  is 
the  specified  "cq.  The  proper  root  x0  is  selected  judiciously  hy 
taking  that  root  for  which  equation  (45)  is  in  closest  agreement  with 
the  approximation  equation  (30) .  The  pressure  pi  is  calculated  hy  use 
of  xc  in  equations  (46)  and  (49) ,  and  is  obtained  from  equation  (47). 

Then  equation  (45)  is  integrated  to  give  the  equation  of  the  line  oh 
(sketch  (d)).  The  intersection  of  that  line  with  that  of  equation  (32) 
locates  point  h . 

In  a  similar  manner,  equations  (45)  through  (50)  are  applied  A 

successively  to  the  ij  states  12  and  23  to  give  values  for  P2  and  c^,  5 

locate  point  c,  and  give  an  expression  for  the  line  cd.  The  inter-  7 

section  of  that  line  and  the  trajectory  of  the  payload  (resulting  from  4 

the  integration  of  eqs.  (20)  and  ( 13) )  gives  point  d  on  sketch  (d) . 

Thus  the  information  needed  to  initiate  the  integration  of  the 
equations  of  motion  and  that  needed  to  estimate  the  regime  of  validity 
of  the  simple  wave  model  has  heen  formulated  hy  hoth  a  series  approxima¬ 
tion  and  the  exact  shock  relations.  Finally,  it  is  instructive  to  * 

consider  briefly  the  simplification  that  arises  when  the  strong  shock 
approximation  is  made. 

The  strong  shock  approximation.-  If  the  value  of  c*!  is  small 
enough,  the  quantity  Cci/^i j)  is  small  and  (ref.  12,  p.  1002)  the 
following  relationships  are  valid 


±  Wjj  _  7  +  1  uj  -  uj 
ci  2  ci 


(51) 


Cj\2  =  7(7  -  1)  /UJ  -  Uj\2 

ci J  "2  V  ci  y 


(52) 


and  as  before 


El  =  7(7  +  1)  /5j  -  uA2 
Pi  2  \  ci  J 


shock  ij 


=  ui  ±  Wi(j 


(53) 

(54) 


* 


2V 
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The  application  of  these  relationships  to  obtain  the  approximate  location 
of  points  b,  c,  and  d  in  sketch  (d)  proceeds  the  same  way  as  for  the 
exact  shock  relations  except_that  now  it  is  not  necessary  to  solve  the 
fourth  degree  equation  for  Wij/ci  which  can  be  computed  directly  from 
equation  (5l) •  For  example,  point  b  is  obtained  by  integrating  equa¬ 
tion  (54)  using  equation  (51)  and  combining  the  result  with  equation  (32) 
which  yields  the  simple  expression 


-  _  7  -  1 
yb  "  7  +  1 

r _ 

Xb  y  +  l 


) 


J 


(55) 


Similarly,  the  expressions  for  p2  and  C2  are  derived  very  simply  from 
equations  (53)  and.  (52)  to  be 


“  =  7(7  +  1) 

Pa  “  2c02(7  -  1) 

-  _  7(7  -  1) 
c2 - g - 


if  7  is  1.41 


(56) 


This  relationship  for  p2  is  quite  good  if  cQ  <  1  and  could  be  used 
where  the  series  expansion  (44)  becomes  invalid. 


Auxiliary  Physical  Relationships 


» 


It  is  useful  to  bring  other  physical  considerations  to  bear  on  the 
problem.  The  problem  of  the  motion  of  the  payload  during  impact  can  be 
solved  if  7,  'c0,  K,  m^/mp,  and  pa  are  specified.  The  selection  of 
these  quantities,  however,  may  be  influenced  by  consideration  of  the 
maximum  allowable  number  of  g' s  acceleration,  allowable  stress  in  the 
bag  skin,  initial  inflating  gas  tenperature,  and  other  conditions.  Thus 
the  auxiliary  physical  quantities  should  be  related  to  the  quantities 
needed  for  solving  the  payload  motion. 


The  maximum  acceleration  of  the  payload  up  to  the  end  of  the  simple 
wave  regime  occurs  at  state  2  in  sketch  (d)  and  is,  from  equations  (12) 
and  (8) , 


P2  "  Pa 
1 


(57) 


18 


♦ 


The  circumferential  stress  in  the  hag  skin  is  a  maximum  at  state  3 
because  of  p3.  However,  the  hag  need  not  withstand  this  pressure  as 
long  as  failure  caused  by  it  does  not  propagate  faster  than  the  shock  cd 
(a  nontearing  hag  material)  .  Then  failure  at  pressures  behind  cd  does 
not  affect  the  motion  of  the  payload  in  the  simple  wave  region  and  the 
maximum  pressure  the  hag  must  withstand  is  p2  •  The  corresponding 
circumferential  force  per  unit  axial  length  (for  a  circular  hag)  is 

(Nc)fc  =  r(p2  -  pa)  »  P0r(p2  -  pa)  (58) 

The  maximum  axial  force  per  unit  circumferential  length  occurs  at  the 
top  of  the  hag  at  state  2  and  is  (from  eq.  (k)) 

W,  "  #  (M\  ™ 

where  (du/dt)^  is  given  in  equation  (57)  •  Combining  (58),  (59),  and  (7) 
yields 

“P  («)  =  ‘  ^b  (6°) 

which  shows  that  Nc  will  he  more  than  double  Na  for  a  positive 
acceleration  at  b.  Therefore,  the  maximum  allowable  stress  is  given 
by  equation  (58) 


Nmax  =  P0r(P2  "  Pa)  (6l) 

Application  of  equation  (57)  at  state  2  using  equations  (ll)  and 
yields 

mp  mpnge 

Because 

mb  =  2jtry0ps  (63) 

(where  ps  is  the  mass  of  skin  per  unit  area) ,  equation  (62)  becomes 

2Jb  =  1 

“p  WPs  - 
2y0nge  '  ^ 


(61) 

(62) 


(6k) 
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Combining  equation  (64)  with  equation  (57)  yields  another  useful 
relationship 

2  (^h/^p)  (  ^pia.v/ Ph) 

^  2K(p2  -  Fa) 


If  the  total  mass  of  the  system  is 

m  =  mg+mp+mb  +  bottom  (66) 

(where  %0ttom  is  the  111888  of  the  "bag  "bottom)  and  the  same  material  is 
used  for  the  bag  sides  and  bottom  so  that 


mhr 

=  =  _2_ 


bottom  =  Jtr  Ps  =  gy 


then 


m  L  mP  mP  V  2yoyJ 


Briefly,  to  solve_a  practical  intact  problem,  one  may  proceed  as 
follows:  Specify  7,  pa,  T0,  and  Uq.  Then  calculate  "cQ  from 


Also  calculate  p2,  c2  and  yb  from  equations  (44) ,  (29) ,  and  (33) 

(if  c0  ;>  l)  or  by  the  use  of  equations  (46),  (47) ,  and  (45)  with 
associated  relationships  (if  “q  <  l) ;  or  use  curves  which  will  be  pre¬ 
sented  subsequently  to  determine  p2,  c2,  and  y^.  Specify  Nmax/Pg* 

Use  equation  (65)  to  obtain  a  suitable  set  of  values  of  m^/mp  311(1  K  "by 
trial  and  error,  so  that  for  the  ”0  calculated,  m^/nip  and  K  fall  within 
the  regime  of  validity  of  the  theory.  (Use  of  a  subsequent  figure  of 
regimes  of  validity  will  be  helpful.)  This  guarantees  that  N^y  will 
not  be  exceeded.  Calculate  mg/mp  from  equation  (ll) .  Specify  r/yD 
(with  buckling  considered  as  noted  earlier) .  Calculate  mp/m  from 
equation  (68)  .  Finally,  specify  n  and  calculate  yc  from  equation  (64) 
and  thus  r  is  determined.  The  motion  of  the  payload  is,  of  course, 
determined  when  the  necessary  quantities  are  used  as  inputs  to  solve 
equations  (20)  and  (13)  simultaneously. 
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I  RESULTS 

i  * 

|  The  differentia],  equations  (20)  and  (13)  subject  to  boundary 

I  conditions  (2l)  were  integrated  numerically  using  the  Adams-Moulton 

f  (ref.  13,  p.  200,  eqs.  6.6.2)  predictor  corrector  method  on  an  IBM  704- 

I  electronic  data  processing  machine.  Initial  conditions  for  the  Integra - 

I  tion  were  obtained  by  both  the  approximate  and  the  exact  methods  described. 

|  These  are  presented  and  compared  in  the  first  three  figures.  Results 

I  obtained  by  the  actual  numerical  integration  are  presented  in  the  last 

|  three  figures. 


Figure  1  is  a  comparison  of  the  first  two  methods  of  computing  p2>  A 
It  is  seen  that  for  large  values  of  c^,,  either  the  exact  shock  relations  5 
or  the  series  expansion  method  (eq.  (44-))  gives  the  same  value  of  p2.  7 
At  c0  =  1,  the  series  method  gives  a  value  of  p2  approximately  10  per-  4 
cent  lower  than  the  exact  method.  For  this  reason,  the  series  method 


should  not  be  used  for  values  of  c‘0  <  1.  However,  for  values  of  "cq  <  1, 
the  strong  shock  approximation  (eq.  (56))  is  sufficiently  accurate  that 
it  cannot  be  distinguished  from  the  exact  shock  relations  on  the  figure. 
Since  p2  was  computed  assuming  that  u2  =  -1  (i.e.,  line  ab  in 
sketch  (d)  is  straight),  the  criteria  derived  for  the  validity  of  the 
wave  model  apply  to  the  figure.  Thus  the  regime  for  p2  <  10  ("Cq  >  l)  is 
not  valid  if  pa  =  0.  For  p2  >  10  Cc0  <  l)  the  wave  model  is  valid 
for  pa  >  0. 

In  many  places  in  the  analysis,  the  quantity  yfc  appears  in  physical 
relationships.  Like  p2,  it  has  been  computed  by  the  series  expansion 
equation  (33  )>  exact  shock  relations,  and  the  strong  shock  approximation 
equation  (55)  and  is  shown  in  figure  2.  The  restrictions  on  the  applica¬ 
bility  of  the  curves  to  the  left  and  right  of  ~cQ  =  1  noted  here  are  the 
same  as  those  on  the  p2  curves  of  figure  1.  In  its  range  of  applica¬ 
bility,  the  series  expansion  yields  good  agreement  with  the  exact  shock 
relations  for  y^,  and  the  two  methods  yield  almost  Identical  results  for 
values  of  c'0  >  4,  The  strong  shock  approximation,  however,  gives  a 
value  of  y^  independent  of  Cq  as  shown. 

I  The  quantity  "c2  required  for  the  integration  of  equation  (20)  is 

I  presented  in  figure  3*  Like  p2  and  y^,  it  has  been  computed  by  the 

|  exact  method,  the  series  expansion  equation  (29),  and  the  strong  shock 

|  approximation  (56),  and  has  the  same  limitations  imposed  above.  The 

series  and  the  exact  methods  give  very  good  agreement,  but  the  strong 
shock  approximation  gives  a  crude  single  value  of  c"2. 

The  regime  of  validity  of  the  simple  wave  model  for  pa  =  1  is  shown 

!in  figure  4.  The  lower  limit  for  the  series  solution  is  at  cQ  ■  1  and 
corresponds  to  the  limit  on  use  of  series  shock  relations.  Thus  the 
series  shock  relations  are  valid  in  most  of  the  regime  of  validity  of  the 
simple  wave  model.  The  upper  limits  on  "S 0  (mentioned  earlier)  for  the 
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various  values  of  m^/mp  correspond  to  the  situation  in  which  the  second 
»  shock,  cd  in  sketch  (d),  just  gets  to  tne  payload  as  it  comes  to  rest. 

It  is  an  approximate  limit  on  the  regime  of  the  simple  wave  model. 

A  corresponding  plot  of  the  regime  of  validity  for  the  case  pa  «*  0 
is  shown  in  figure  5.  The  same  comments  apply  except  that  part  of  the 
upper  limit  is  determined  hy  the  value  c0  «  1  to  insure  that  p2  ^  10 
in  accord  with  the  criteria  developed  in  the  analysis.  Thus,  application 
of  the  present  wave  model  to  impact  without  atmosphere  is  limited 

essentially  to  situations  in  which  0  <  cl  <  1. 

—  u 

The  solutions  of  the  differential  equations  (20)  and  (13)  are 
presented  in  figure  6  for  one  impact  condition.  The  solutions  were 
obtained  by  use  of  both  series  and  exact  methods  for  a  value  of 
Co  =  1.87b  (which  is  not  expected  to  give  outstanding  agreement).  The 
differences  in  y^  and  ym^n  between  the  two  methods  are  about  2  and 

5  percent,  respectively.  The  difference  in  p2  is  about  3  percent. 

The  velocity  computed  by  the  series  method  decays  at  about  the  same  rate 
as  that  computed  by  the  exact  method  but  lags  by  about  5  percent  of  the 
initial  speed  at  any  given  instant.  Thus  the  agreement  between  the  two 
methods  is  satisfactory. 


CONCLUDING  REMARKS 


A  theory  of  impact  cushioning  by  use  of  cylindrical  gas  bags  has 
been  developed.  The  theory  is  applicable  to  high-speed  impact  conditions 
because  unsteady  gas  dynamic  phenomena  have  been  incorporated  in  the 
analysis  by  the  use  of  a"  simplified  mathematical  model  for  the  waves  in 
the  gas.  The  influence  of  the  bag  skin  forces  on  the  motion  of  the  pay- 
load  has  also  been  accounted  for  in  the  analysis.  The  shock  waves  in 
the  gas  are  treated  by  exact  relationships,  the  strong  shock  approxima¬ 
tion,  and  a  series  expansion.  The  last  two  treatments  lead  to  relatively 
simple  useful  physical  relationships  for  analyzing  the  impact  motion. 

The  methods  are  compared,  and  the  regime  of  validity  of  each,  as  well 
as  the  regime  of  validity  of  the  basic  wave  model,  is  evaluated.  Curves 
of  some  physical  results  are  presented  which  can  be  used  with  equations 
presented  for  calculation  purposes. 
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Figure  2.-  Height  above  ground  where  shock  first  strikes  payload;  7 


Figure  4.-  Regime  of  validity  of  simple  wave  solution;  pa  =  1,  7  =  1.4l. 


Figure  6.-  Comparison  between  exact  and  series  methods;  7  =>  1.41, 
c0  =  1.874,  K  =*  2,  ffib/mp  »  1,  pa  =>  1. 
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